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, In the present paper we study the role of particle interactions on the evolution of a high energy 

• beam. The interparticle forces taken into account are due to space charge alone. We derive the 

collision integral for a charged particle beam in the form of Balescu-Lenard and Landau and con- 
sider its further simplifications. Finally, the transition to the generalized kinetic equation has been 
accomplished by using the method of adiabatic elimination of fast variables. 



I. INTRODUCTION. 



In most of the works so far, dedicated to the study of beam plasma properties the effect of interparticle collisions 
has been neglected. In many important cases this is a sensible approximation giving satisfactory results, yet one has 
to elucidate the limits of validity of " collisionless beam" approach and to investigate the role of collision phenomena 
in beam physics. Collisions are expected to bring about effects such as thermalization, resistivity, diffusion etc. 
that influence the long term behaviour of charged particle beams. The reasoning commonly adopted for employing 
the "collisionless beam" approach is that characteristic beam-plasma frequencies are much greater than collision 
fS frequencies for a large number of situations in beam physics. Such an assumption is not based on stable physical 
grounds as pointed in Q]. 

The term " collisionless beam" means that interactions between particles giving rise to dissipation and hence leading 
to establishment of equilibrium state are not taken into account. In a number of cases involving reasonable approx- 
■^j- ■ imations it is sufficient to compute the macroscopic characteristics (charge and current densities) in relatively big 
volume elements containing a large number of particles. As a result interaction manifests itself in the form of a 
mean, self-consistent field thus preserving the reversible character of the dynamics involved, and leading to the time 
On ■ reversible Vlasov's equation. 

The notion of " collisional beam" usually conceived as the counterpart of " collisionless beam" implies that dissipation 
due to redistribution of beam particles is taken into account, resulting in additional term (in the form of Landau or 
Balescu-Lenard) in the kinetic equation. In a sense, Landau and Vlasov approximations correspond to two limit cases: 
namely the Landau collision integral takes into account interactions that determine dissipation while the effect of the 
• i— i ' mean, self-consistent field is not included into the physical picture involved. On the contrary, the latter is the only 
way interactions manifest themselves in the Vlasov equation, leaving however the question about the role of collisions 
near particle-wave resonances unanswered. The Balescu-Lenard approximation lies somewhat in between Landau and 
O ,. Vlasov limit cases with the due account of dynamic polarization of the beam, that is a more complete inclusion of 
ILJ ■ collective effects resulting from interactions between charged particles. 

In the present paper we derive the collision integrals for charged particle beams. The transition to the unified kinetic, 
hydrodynamic and diffusion description of particle beam propagation embedded in the generalized kinetic equation 
P] is further accomplished building on the concept of a coarse-grained hydrodynamic picture. The latter implies the 
existence (and their proper definition) of characteristic spacial and temporal scales typical for the hydrodynamic level 
of description Q , || . Within the elementary cell of continuous medium thus defined it is naturally assumed that local 
equilibrium state is reached. This state is further described (defining the drift and diffusion coefficients in coordinate 
space) by the method of adiabatic elimination of fast variables, widely used to match the transition to Smoluchowski 
equation jij . The granulation of phase space with the due account of concrete structure of continuous medium results 
in additional collision integral in the kinetic equation, thus describing the dissipation caused by spacial diffusion of 
the distribution function and redistribution of particle coordinates |§ , ■ 

The generalized kinetic equation makes it possible to build an unified picture of non equilibrium processes on kinetic 
and hydrodynamic scales without involving a perturbation expansion in Knudsen number M. It can be shown that 
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the set of hydrodynamic equations for cold beams put in appropriate form is equivalent to mesoscopic quantum-like 
description of particle beam propagation |J], §. 

The scope of the presentation given in the paper is as follows. In Sections II and III we formulate and solve the 
equation for the fluctuations of the microscopic phase space density in the case of a space charge dominated high 
energy beam. The solution obtained provides the grounds to find explicitly the collision integral in the form of 
Balescu-Lenard. Sections IV - VI deal with the various forms and simplifications of the collision integral. In Section 
VII we derive the additional Fokker-Planck term in the generalized kinetic equation. Finally, Section VIII presents 
the conclusions of our study. 

II. AVERAGED MICROSCOPIC EQUATIONS. 

In a previous paper ^ we derived the equation for the microscopic phase space density with a small source, taking 
into account the proper physical definition of continuous medium. It was the starting point in the transition to the 
generalized kinetic equation for the one-particle distribution function. The equation for the microscopic phase space 
density reads as 
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where N = N (x, p'*^; 0) is the true microscopic phase space density written in the variables 



x = (x, z, a) 



(2.2) 



They are related to the canonical coordinates x = (x, z,cr) and canonical momenta p = (p x ,p Zl h) through the 
following equations 
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where u — (x, z) and all other notations are the same as in Ref. 0]. In particular the following designations 
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e s = (0,0,1) (2.4e) 
have been introduced in equation (2.1), while N (Si,p^;8) is the smoothed microscopic phase space density 
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with a smoothing function (x| p). 

The next step consists in averaging the Klimontovich equation (2.1) over the relevant Gibbs ensemble with using 
the definition of one-particle distribution function tt 



iv(x,p( fe );0)) =nf (x,pW;fi' 
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where iV p is the total number of particles in the beam and V is the volume occupied by the beam. By taking into 
account the representation of the microscopic phase space density and the microscopic force in terms of mean and 
fluctuating part 



N = nf + SN ; F (M) = (Fj + SF ; ^7VF (M) j = nf (Fj + (SNSFj 



we obtain the generalized kinetic equation 
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are the collision integrals. It was previously shown ||, j| that the additional collision integral ^(x, p^';S) can 
be cast into a Fokker-Planck "collision term" , where the Fokker-Planck operator acts in coordinate space only. The 
equation for the fluctuating part SN reads as 
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Averaging the Maxwell-Lorentz equations we get 
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V r • (B) = 



V r -<E) = ^p(r;t) ; 
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where 



p(r;t) = ydV fe) /(r,P (fe) ;i) 5 j (r; t) = J dV fc >v/ (r, p«; 
The equations for the fluctuating fields are similar to (2.10) and read as 
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where 



5p(r;t)= J d 3 p^5N(v,p^;t) ; Sj (r; t) - / dV fc >v<WV (r, p«; i) . (2.13) 

Taking divergence of the first of equations (2.12) and utilizing the last one, it can be easily seen that the continuity 
equation for fluctuating quantities holds 
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It should be pointed out that the microscopic electromagnetic fields depend on the coordinates x = (x, z, a) 
through the microscopic phase space density N written in these coordinates. The rest of this section is dedicated to 
the derivation of some useful relations, needed for the subsequent exposition. Consider the simple change of variables 
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Continuing further we use the relations 



and finally get 



As far as 



u = u + rfD u ; j#>=p«+5^i ; a = a+ ]T (p u V u - } . 

u—(x,z) 

h {k) = s{k) + J_ 



i I «2;>(fc) 

d 3 rd 3 p (fe) = p 3 (1 + x . K) 1- Po 'I dxdzdad^ dftV d rfk) _ 



S S3 y / l + 2?y( fe ) + /?2^(fc)2 ~ x 



for ft « 1 we obtain 



d 3 rd 3 p (fe) = p 3 (1 + x . K) d 3 2d 3 p (fc)_ ( 215 ) 

Thus, integration in the expressions for the charge and current density 

Sp(r;t) = j d 3 p^8N (x,p^;8^) ; <5j (r; t) = ^ d 3 P < fe KSN (x, p W ; 6>) . (2.16) 
goes approximately over the new kinetic momenta p( fe ). 

III. SPECTRAL DENSITIES OF FLUCTUATIONS. 

In order to determine the collision integral (2.8) we have to solve equation (2.9) governing the evolution of fluctua- 
tions SN. Under the assumption that fluctuations are small the second term on the right hand side of equation (2.9) 
can be neglected 



SN (x; 6>) = 



= -ni2V p -(<SF(x;0)/(x;0)). (3.1) 

The small source in the initial equation (2.9) has been dropped off as non relevant for the dynamics of small-scale 
fluctuations. The term containing the mean force in equation (3.1) can be neglected. This is justified when calculating 
the small-scale fluctuations if 

( 2 Q 2n 2 £ok B T\ 

v *»<*>° \^ = 7^r ; rD = ^r)- (3 - 2) 

Here T is the temperature of the beam, uj is the angular frequency of synchronous particle, v x ,z,cr stands for the 
betatron tunes in the two transverse planes as well as for the synchrotron tune. Furthermore u p is the beam plasma 
frequency and rjj - the Debye radius. It is worthwhile to note that the physical meaning of Debye radius for particle 
beams is somewhat different from that commonly used in plasma physics. In fact Debye radius is an equilibrium 
characteristic of the beam, indicating the exponential decay of the self-field, needed to self-maintain this equilibrium 
state. 

The contribution of small-scale fluctuations can be better extracted if a small source proportional to A is introduced 
into the left hand side of (3.1) 
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In going over from equation (3.1) to (3.3) the left hand side has been represented in terms of the variables x = (a;, z, a) 
The general solution of the above equation can be written as 



SN (x, p (fe) ; fl) = SN S (x, P (fc) ; 0) + SN md (x, P (fe) ; 0) 



(3.4) 



where 8N md is a generic solution of (3.3), while SN S accounts for the discrete structure of the beam as a collection 
of particles. The latter can be determined from 
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with the initial condition 



(SN S (X;6)SN S (X x ;0)) = nS (x - x x ) S fpW - p^ fc) ) / (x, p« ; (9 



(3.6) 

When small-scale fluctuations are computed / (x, p( fc ) ; 0) can be considered a smooth enough function (not varying 
considerably) and (5N S (X; 8) 5N S (Xx; depends on 9 — Q\ and x — xi only. Introducing the Fourier transform: 



(SN (X; 6) SN (X x = (SNSN) 
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we cast equation (3.5) into the form 
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Further we introduce the one-sided Fourier transform in the time domain 
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Multiplication of equation (3.5a) by e iWT and subsequent integration on r yields: 
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Using the equation 
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relating the one-sided and two-sided Fourier transform we get 
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The definition of Dirac's 5-function 
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gives finally 
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To obtain an arbitrary solution of equation (3.3) we perform the Fourier transform 
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What remains now is to compute the spectral density of fluctuating force SF. In doing so we consider an arbitrary 
function F (x; 9). Let the same function, written in the variables r = (x, z, s = R9) and t be F r (r; t). Further we have 



F r (r;t) = — ^— j / W 3 mF r (f;m)e 
(2tt) 
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where the following change of variables 



j/ = v n k n 



has been introduced. Therefore the relation we are looking for reads as 
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Fourier analysing equations (2.12) we find 
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Let us represent the electromagnetic fields as a sum of longitudinal and transversal components 



im-SE r = — Sp r . 



SE r = 5E| + SE^r (m x <5Ej! = ; m-SE^r = o) 
and further simplify the problem by considering 



From the continuity equation (2.14) we get 



Sjr = v e s 5p r . 
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v 



and using (3.13) and (3.16) we conclude that m = k. Thus we obtain 
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Retaining leading terms only, we write the fluctuating force 5F as 



SF (to, k) = [<5E + v (e s x 5B) 
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Integrating equation (3.11) on p( fc ) we obtain 



zg 2 k 



e /%tfE o (P - dlkl) 



Sp (w,k) . 



.SpU-.k) = U:k) - inll I rf 3 p( fc ) - V ^ k ( ~^ A -£F>,k) 



and eliminating <5F (w, k) with (3.19) in hand we get finally 
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where 



ZaPlllE (fc 2 - /3 2 fc 2 ) J u>- Kk-v + tA 

is the dielectric susceptibility of the beam. Thus for the spectral density of the fluctuating force we have the following 
expression: 

IV. COLLISION INTEGRAL IN THE FORM OF BALESCU-LENARD. 

According to (2.8) the collision integral is given by 

J col (x, P « -e)=~%- (sFSn) (x, P « , 9; x, p « , d) . (4.1) 

We shall express the right hand side of (4.1) in terms of the spectral densities of fluctuations SF and SN. Let T (x; 9) 
and Q (xi; #i) be two random functions. The second moment in the variables x — xi, 9 — 9\ can be written as 

{TG) (x,0;xi,0i) - (TQ) (x,0;x-xi,0-0i) = 

= — ^ / dwd 3 k (?<f) (w, k; x,0) cxp {i [k- (x - xi) - u (9 - 6»i)]} . (4.2) 
(27r) J ^ ' 

As far as the second moment is a real function the spectral density obeys 

(JO) (uj, k; x,0) = * (-uj, -k; x,0) . (4.3) 
Letting x = xi, 6 = 6\ in (4.2) with (4.3) in hand wc find 

{TG) (x,0; x, 0) = — ^— j / rfwd 3 k i?e (^) (w, k; x,0) . (4.4) 
(27r) J ^ ' 

Using (4.4) and taking into account only leading terms in SF we rewrite (4.1) as 

J col (x,p {k) ;0)=-*%-J ^Re(s¥sN^ (,,k;x,pW,0). (4.5) 



Utilizing the expressions (3.11) and (3.22) we obtain 



where 
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/ d 3 p«/(x;0)2 7 r ( 5( W -i?k.v). 



(4.6) 



In formula (4.5) representing the collision integral only the real part of ^<$F£/V^ (w,k;x, p( fe ),#) enters. Therefore 
the expression to be substituted back into (4.5) reads as 

Re (dFSN^j (u), k; x, pW , 6») = -nnRS (w - Rk-v) k ^<5F<5F^ k • V p / (x; 0) - 
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where 
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J d 3 p^S (w - .Rk-v) k • V p / (x; 9) . 



e^ ol lE (fc 2 - 02fc2) 

Finally, the collision integral (4.5) can be written in the form of Balescu-Lenard as 

7 {BL) h nW-fl) - y . f dV ^ k) 5(1, y k-viU 



(4.8) 
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• / (x, pf» ; 0) V p / (x, p« ; o) / (x, p« ; *) V Pl /(x,pf»; *) 
The collision integral (4.5) can be put in an equivalent form of a nonlinear Fokkcr-Planck operator 



(4.9) 



V p / (x, 



A' Bi »/(x,p w ;fl) 



where the drift and diffusion coefficients 



(4.10) 



& BL ^ = ttR 2 { ^-±5 ( W - i?k • v) k f <5F<^ k, 

7 (27r) V /o>,k 



(4.11a) 



e ^E o J (2tt) 4 
depend on the distribution function itself. 



k /me (w, k) 
fc 2 -/? 2 *; 2 |?( w ,k)| 2 



(4.116) 



V. COLLISION INTEGRAL IN THE FORM OF LANDAU. 

The dielectric function (3.21) depends on the distribution function and consequently the corresponding kinetic 
equation with the collision integral in the form of Balescu-Lenard is extremely complicated to solve. Thus one should 
seek reasonable ways for further simplifications. First of all we shall determine the equilibrium state described by 
/o (x, p( fe );#) and satisfying 
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^ + R (v • V a ) + i? (f • V p ) /o (x, pW; fl) = 0. 



(5.1) 



It can be easily checked that equation (5.1) has a solution of the form 



'• I - —(k) a\ c ( ^ x 



(5.2) 



where 



2 ^ ~ Pu 



(it = x, z) , 



(5.3a) 



2J CT 



ffk)2 



+ X \d - a s 



R 



tan<&. 
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In the above expressions a, (3 and 7 are the well-known Twiss parameters 

da u G u dp u d-y u 2G U 

— = — Pu - Rlu ■ -^- = -2Ra u ; — = —a v 



= 1, 



while 



1 A^o K 

2ttR 2 PIE K 



cos 



(j9„ = A" 1 ) ; ^ - i? 2 /C 2 A 2 , 



k - being the harmonic acceleration number, <!> s - the phase of synchronous particle, v s is the synchrotron tune and 
/So- can be interpreted as the "synchrotron beta- function" . The quantities e x , e z and e CT are related to the transverse 
and longitudinal beam size and are referred to as equilibrium beam emittances. To describe a local equilibrium state 
(see next section) one can formally choose the equilibrium beam emittances e x , e z and e a proportional to the beta- 
functions by a universal scaling factor e/R characterizing the equilibrium state. Let us recall that at local equilibrium 
all the parameters of the distribution are allowed to depend on coordinates and time [Q , which is consistent with the 
specific choice above. Further, by specifying the generic function (5.2) for slowly varying beam envelopes we find 



3/2 



exp 



2R(j x f3- 1 + J z f3- 1 + J a l3- 1 ) 



2c 



(5.4) 



The equilibrium beam emittances e x , e z and e a are related to the temperature of the beam through the expression 



k B T 



(5.5) 



In order to obtain the collision integral in the form of Landau we consider e(Rk ■ v,k) = 1 in equation (4.9) and 
simultaneously take into account the effect of polarization by altering the domain of integration on k for small k. 
As far as the large values of k are concerned the upper limit of integration can be obtained from the condition that 
perturbation expansion holds. To proceed further it is convenient to change variables in the Balescu-Lcnard kinetic 
equation according to 



-sign (/C) 



V 



k n 
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This means that the canonical coordinate a has been transformed according to a — ► <7-i/|/C|, and in order to retain 
the hamiltonian structure of the microscopic equations of motion the a - component of the force should also be 
transformed as F a — ► — sign (JC) F a . Taking into account the fact that the Balescu-Lenard collision integral is 
proportional to the square of the fluctuating force we can write 



j{BL) / x 6 ( fe) . A _ nq*nR ~ f dVy^ / (fc) (fc) 
J col ^x,p ,°)- £ 2 ME2 J p J (2 ^3 d(kp k Pl 
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where 



Handling the integral 



/ (x, pf } ; 9) %f (x, p( fc > ■6)-.f (x, p(*> ; 6) V P1 / (x, p[ k) ; O) 



I L (g) = / d 3 k 



kk 



{kl + kl + klhl\K\)- 



r<5(k-g), 



(g = P W -p£ fc) ) 



(5.6) 



(5.7) 



by choosing a reference frame in which the vector g points along the a - axis, and using cylindrical coordinates in this 
frame we find 



2-7T OO 



I L (g) = J dk ± k ± Jd<£ J dk a 8{k a g) 



1 



-oo 



{ki+kihi\K\y 



k± cos <fr \ / k± cos $ 
* | k± sin $ fc_L sin $ 



fe_D 



(5.8) 



As was mentioned above in order to avoid logarithmic divergences at both limits of integration on k± in (5.8) we have 
altered them according to 



kn 



1 



AlTEokBT 
loQ 2 



Thus the Coulomb logarithm C is defined as 



(5.9) 



£ = ln 



47T 



(eok B T) 3/2 



The tensor Ij, (g) can be evaluated in an arbitrary reference frame to give 

Finally the collision integral (5.6) can be represented in the form of Landau as 

q 4 nRC 



3% (x,p<*);(?) = 



/ 



V p - /d 3 pl fe) G L (g) 



(5.10) 



(5.11) 
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/ (x, pf» ; 9) %f (x, p« ; 9) - / (x, p« ; e) V Pl / (x, p< 



(*). 



(5.12) 



where 



G i( g) = i(l-f 
5 V 9 



is the Landau tensor 



VI. THE LOCAL EQUILIBRIUM STATE AND APPROXIMATE COLLISION INTEGRAL. 

The local equilibrium state is defined as a solution to the equation 



J col (x,pW;0) =0, 



(5.13) 



(6.1) 



where the collision integral is taken either in Balescu-Lenard or Landau form. This solution is well-known to be the 
Maxwellian distribution 



•5 0*" - " 



(6.2a) 



/" d 3 xp(x;0) = F, 



(6.26) 



where p (x; 0), e (x; 0) and u (x; 9) arc functions of x and 8. It should be clear that the local equilibrium state is not 
a true thermodynamic equilibrium state, since the latter must be homogeneous and stationary. To prove that the 
distribution (6.2) is a solution of (6.1) when the collision integral is taken in Landau form (5.12) it is sufficient to take 
into account the obvious identity 

G L (a) • a = a T • G L (a) = 0. (6.3) 
Next we note that the Landau collision integral (5.12) can be written as a nonlinear Fokker-Planck operator 



where 



J%> (x, P W ; 0) = B [V, • (V ■ V p ) - V p • A] / (x, p W ; 9 



(6.4) 



B = 



C 



(6.5a) 



V= I d*p[ k) G L ( g )f(p[ k) ) 



A = J d 3 p (k) G L (g)-%j(p[ k) 



(6.56) 



Our goal in what follows will be to match the transition to the unified kinetic equation. To approach this it is sufficient 
to compute the drift and diffusion coefficients (6.5b) using the local equilibrium distribution (6.2). A more systematic 
approximation methods using the linearized Landau collision integral can be found in Q. Going over to the new 
variable 
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we write 



*(C-|/ 



d 3 dG L (g) /, (CO = A (C) G L (6p) + A 1 (C) 



<5p<5p 



where g = C — Cj and _4 , -4i are functions of the modulus of the vector C 



(6.7) 



A 1 (C)=Sp-V-6p=^ I f / :i C L rfp-G L (g) <*p/, ; (C,). 



(6.8a) 



^0(^ = 1^(^-^(0, 



(6.86) 



A(Ci) 

ic-cj- 



(6.8c) 



To compute the integrals (6.8a) and (6.8c) we use spherical coordinates in a reference frame in which vector C points 
along the a - axis. We find 



2tt 1 



^T 2 V :§ / dcriC i f M f dcos0 ~ 



-ci 



2tt 1 



Ai (C) = y ddCf j d$ J rfcosG^ 



2 (c 2 -cCicose) z 



c 2 - 



where we have used g • C = C 2 — CCi cos 9. Changing variables in the above integrals according to 



. 9 2 = C 2 + C 1 2 -2CCiCose 



dcos0 = 



CCi 



dg 



yields the result: 



s P (v) 



Ap / R 



c+d 



C V 27T6 



dg 



|C-d| 



8p R 
~C V 2^ 



2p /i? 
C7 



erf (C) . 



and similarly 



C+Ci 



Off 



|C-d| 



2 (g 2 -C 2 1+ C 2 



2 1 



c 2 - 
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VII. THE GENERALIZED KINETIC EQUATION. 



The transition to local equilibrium, that is the kinetic stage of relaxation, is described by the Balescu-Lenard or 
the Landau kinetic equation. The latter with due account of the approximate collision integral (6.15) can be written 



f + H(p W -V,)/^(F.V P )/ = J- 



5 vj + v, 



p« - u 



(7.1) 



F = F + (F) 

It is well-known 0] that the kinetic equation (7.1) is equivalent to the system of Langevin equations: 



rfx 



dpW 1 



;(fe) 



u) +RF 



d6 dO @rel ^ ' V R$rel 

where £ (9) is a white-noise random variable with formal correlation properties 



(7.2) 







(M0Kn(0i)) = 2<W(0-0i). 



(7.3) 



The generalized kinetic equation (2.7) describes the evolution of the beam for time scales greater than the relaxation 
time 9 re i. In order to determine the additional collision integral J (x, p^;6j we use the method of adiabatic elim- 
ination of fast variables, which in our case are the kinetic momenta p( k \ In the limit of small times 9 re i (compared 
to the time scale of physical interest) the second equation (7.2) relaxes sufficiently fast to the quasi-stationary (local 
equilibrium) state for which dp^ /dO — ► 0. Thus we find 



;(*0 



u + R0 re iF + 



~~R~ 



(7.4) 



and substituting this into the first of equations (7.2) we arrive at 



dx 

~d~e 



Ru + R 2 6 rel F + y/eR9 rel £(6) . 



(7.5) 



The above equation (7.5) governs the evolution of particles within the elementary cell of continuous medium, where 
local equilibrium state is established. Such a coarse-graining procedure gives rise to the additional collision integral in 
the generalized kinetic equation (2.7). The latter follows straighforwardly from (7.5) and can be written in the form: 



j(%pW;0) =R0 re i{Vx- [e(x;0)V*] - R (v* • f) } / (x,p<*>; 9 



(7.6) 



VIII. CONCLUDING REMARKS. 



In the present paper we have studied the role of electromagnetic interactions between particles on the evolution of 
a high energy beam. The interparticle forces we have considered here are due to space charge alone. Starting with 
the reversible dynamics of individual particles and applying a smoothing procedure over the physically infinitesimal 
spacial scales, we have derived a generalized kinetic equation for kinetic, hydrodynamic and diffusion processes. 

We would like to point out an important feature of the approach presented in this work. The irreversibility of beam 
evolution is introduced at the very beginning in the initial equation (2.1) for the microscopic phase space density. 
Smoothing destroys information about the motion of individual particles within the unit cell of continuous medium, 
hence the reversible description becomes no longer feasible. Details of particle dynamics become lost and motion 
smears out due to dynamic instability, and to the resulting mixing of trajectories in phase space. 
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The collision integral for a high energy beam has been derived (Sections IV and V) in the form of Balescu-Lenard 
and Landau. This collision term scales as E~ 6 (E a is the energy of the synchronous particle) which comprises a 
negligibly weak dissipative mechanism for high energy beams. 

To accomplish the transition to the generalized kinetic equation the Landau collision term has been simplified by 
linearizing it around the local equilibrium distribution. The latter suggests a close relation between equilibrium beam 
omittance and the temperature of the beam. 

Finally in Section VII we have derived the additional dissipative term due to the redistribution of particle coordi- 
nates. This has been achieved by applying the method of adiabatic elimination of fast variables (kinetic momenta). 
The physical grounds for this application is provided the fact that within the physically infinitesimal confinement 
the relatively slow process of smear in configuration space is induced by the sufficiently fast relaxation of particle 
velocities towards a local equilibrium state. It maybe worthwhile to note that a more systematic approach involving 
the projection operator technique Q could be used to derive the additional collision integral in the generalized kinetic 
equation. 
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